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ABSTRACT 

We analyze the gamma-ray emission from 9 high latitude, translucent molec¬ 
ular clouds taken with the Fermi Large Area Telescope (LAT) between 250 MeV 
and 10 GeV. Observations of gamma-rays allow us to probe the density and 
spectrum of cosmic rays in the solar neighborhood. The clouds studied he within 
~270 pc from the Sun and are selected from the Planck all-sky CO map. Gamma- 
rays in this energy range mostly result from cosmic ray interactions with the 
interstellar medium, which is traced with three components: HI, CO, and dark 
gas. Every cloud is detected and shows signihcant, extended gamma-ray emission 
from molecular gas. The gamma-ray emission is dominated by the CO-emitting 
gas in some clouds, but by the CO-dark gas in others. The average emissivity 
and gamma-ray power law index from HI above 1 GeV shows no evidence of 
a systematic variation. The CO-to-H 2 conversion factor shows no variation be¬ 
tween clouds over this small spatial range, but shows signihcant variations within 
each cloud. The average CO-to-H 2 conversion factor suggests that the CO-dark 
gas is molecular as opposed to optically thick HI. 

Subject headings: cosmic rays — gamma-rays: ISM — ISM: clouds 


1. Introduction 


Molecular clouds at high galact k latitude represent an abundant source of gamma-rays 


due to interactions with cosmic rays flDigel et ah lll996l:lAbdo et ah ll2010l:lAckermann et al. 
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2012a ^0; Marti et al. 2013 '). High latitude clouds are primarily low mass and harbor 


little to no active star formation (see iMcGehee l2008l for a review). They are also often 


relatively isolated from localized cosmic ray acceleration sites such as supernova remnants 
or OB associations. Thus, their gamma-ray emission should reflect the steady-state cosmic 
ray density and spectrum of the surrounding region in the Galaxy. These clouds are all 
nearby with most hay ing a distance d < 350 pc for |6| > 25° given a scale height of 150 pc 
flMagnani et al. Ill996l ). High latitude clouds therefore represent potentially pristine probes 
of the cosmic ray spectrum in the solar neighborhood. 


Diffuse gamma-ray emission comes from a combination of cosmic ray interactions with 
the interstellar medium (ISM), inverse Compton scattering of ambient radiation by cosmic 
ray electrons and positrons, and extragalactic diffuse emission. Cosmic ray interactions with 
the ISM produce primarily GeV gamma-rays via proton-proton collisions which lead to the 
production of neutral and charged pions. Neutral pions decay directly into gamma-rays while 
the charged pion species ultimately decay into electrons and positrons. Cosmic ray leptons 
interact with gas to create gamma-rays mainly via bremsstrahlung emission. The ISM is 
effectively transparent to cosmic rays, so the gamma-ray emission is sensitive to the total 
gas column density regardless of dust properties or gas state (HI, HII, H 2 ). 


The necessary likelihood modeling of Fermi Large Area Telescope (LAT) data incor¬ 
porates contributions to the gamma-ray emission from gas traced b y 21 cm HI emissio n 


and the integrated 2.6 mm CO(J = 1 —)■ 0) line emission (ITco) flLebrun et al. 1198311 . 


Hco is assumed to be directly proportional to the column density of molecular hydro- 
gen: A^fH^) = Xr.oW r.n ■ In fact, gamma-rays have been used to calibrate Xco (e.g.. 


Bloemen et al. 119841 1. However, likelihood models of gamma-ray emission based solely 


on the distribution of these two species exhi bit significant residual gamma-ray emission 
flGrenier et al. 1120051: lAckermann et ah 11201 Tail . This excess gamma-ray emission is assumed 
to trace unseen molecular or atomic gas, which has been called dark gas. In cold molecular 
clouds, the dark gas mass can be a signih cant fraction of the total gas mass flAde et al. 112011 


Paradis et al. 

2012; 

Pineda et al. 

2013) 


Dark gas is expected in photodissociation regions and diffuse clouds exposed to the inter¬ 
stellar radiation held. Because GO self-shields less efficiently and has a lower dissociation en¬ 
ergy than H 2 , CO dissociation occurs deeper int o a cloud (i.e., to higher Ay) Therefore, be 


tween 1 < Ay < 5, CO fails to trace H 2 linearly (Ivan Dishoeck fc Black 


1988 


Wolhre et al. 


2010 1. Clouds in this range of Ay are classihed as “trans lucent” (Ivan Dishoeck fc Black 


19881 ). and represent most molecular clouds at high latitude flMagnani et al. 


19961). Translu¬ 


cent clouds tend to be smaller and less den se than giant molecular clouds, but should be 


much more numerous flMagnani et al. 11198511 
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A lingering problem with modeling gamma-rays from high latitude clouds had been that 
existing CO maps were unavailable, sparse, or i ncomplete. The Planck CO rn ap is the hrst 
all-sky map of the CO( J = 1^0 ) emission line nPlanck Collaboration I l2014b|l . We use this 
CO map in a complete, flux-limited survey of high latitude molecular clouds with the Fermi 
LAT. In this paper, we describe the methods and systematic uncertainties of the survey and 
report on gamma-ray observations of 9 high latitude, translucent molecular clouds. Every 
cloud has detectable gamma-ray emission consistent with maps of either CO, dark gas, or 
both. We also report the photon index of each cloud, as it reflects the incident cosmic ray 
spectrum. This work presents the initial results from the full survey. 


2. Source Selection 


For this pilot survey, we identify well-known molecular clouds with |6| > 25° from 
previous surveys to analyze in gamma-rays. Clouds at such latitudes are all nearby, which 
ensures the gamma-ray flux is high enough to study. Large clouds at low er galactic latitude. 


such a s Orion, have already been extensively studied with the Fermi LAT flAckermann et aL 

201230 ) • 


Torres et al. (12005) presented a study on the possible garnma-ray emission based on a 


number of CO surv eys flMagnam et al. Ill996l: iHaTtmann et al. Ill998l: iMagnani et al. 112000 


Dame et al. 120011) . We find candidate clouds from these surveys and identify them in 


the Planck CO map. We choose several of the brightest high latitude clouds from these 
surveys. Three other, fainter clouds are chosen to explore the low ILco limits for gamma- 
ray detection. Two bright clouds, G313.1-28.6 and G315.1-29.0, are identified via visual 
inspection of the Planck GO map. These were observed in an earlier GO catalog to wards 


( 120011 ) 


namea v-nair 

Torres et al. 

laeieo 

(2005 

n-gjasL 11 liviizuno 

), nor mapped by 1 

zuuii. Dur 

Dame et al. 


In addition, we choose a region devoid of large-scale g as and dust emission by visua l 
inspection of the Planck GO map and the color excess map oi lPlanck Gollaboration fj2014all . 
This region, centered on Galactic coordinates (£, 6) = (250°, 30°), is used to test the false 
detection rate of gamma-ray emission from GO or dark gas. 

Some relevant properties of the chosen clouds are listed in Table [U The areas of the 
clouds are calculated from the GO extent in the Planck map. Distances to the clouds 


are taken from the lit erature and were derived from color excess fjSchlaflv et al. I2014J: 


Lallement et al. II2014J ). Many, but not all, of the clouds have masses estimated in prior 





























































-4- 


surveys, (e.g., 


Magnani et al. 11 19961) . 


The locations of the clouds relative to the Solar System can be seen in Figures [H [21 and 
El Most of the clouds lie below the galactic plane, as shown in Figures [2] and El but represent 
a large range of Galactic longitudes seen in Figure [H 


3. Gamma-ray Analysis 


The LAT on the Fermi Gamma-ray Space Telescope is a pair-tracking telescope, sensitive 
to gamma-rays between 20 MeV and 300 GeV. The tracker is surrounded by anti-coincidence 
detectors to distinguish between cosmic ray and gamma-ray events. The photon localization 
strongly depends on the photon energy; at 1 GeV, the 68% containment radiu s is 0?8, 
decreasing with energy to roughly 0?2 at around 10 GeV flAckermann et ah Il2012el) . 


We use data of the entire sky from the Fermi LAT between August 4, 2008 and June 19, 
2013. We use the Fermi science tools (v9r27pl) available from the Fermi Science Support 
Genter 0, utilizing the P7_V6 instrument response function. When selecting the data, we 
consider both front and back converted photons in the “source” class. We select data from a 
10° radius around the chosen coordinates and between 250 MeV and 10 GeV. These energies 
are chosen to maximize both source localization and photon statistics. Including photons 
between 10 and 100 GeV does not improve the signihcance of the detection, as will be 
explained in Section jH We exclude photons with incidence angle > 100° from the zenith 
and any time the spacecraft rocking angle exceeds 52°. These constraints remove most 
gamma-ray contamination coming from the Earth’s limb. 


We perform a binned likehihood analysis, selecting the data which lie inside a 14° x 14° 
square centered on the region of interest (ROI) center. The ba s ic pro cedure for the likeli¬ 
hood analysis of gamma-ray data is described in iMattox et al. (119961) . To evaluate source 
detection and model signihcance, we consider the test statistic (TS), which is proportional 
to the difference of the log likelihoods of two different models: 


TS = -2(^lnCA-\nCB), ( 1 ) 

where Ca and Cb are the likelihoods for two models we are comparing. The TS represents 
the signihcance of model B over model A. 


^FSSC: http://fermi.gsfc.nasa.gov/ssc/ 
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To create the model, point sources are taken from the Fermi 2 year catalog flNolan et ah 


2 OI 2 I) version 6 (2FGL). For those clouds with point sources from the 2FGL coincident with 
the GO emission, we remove the point source from the model. Sources outside of 7° from 
the ROI center have all free parameters hxed for the htting procedure, while closer sources 
only have their spectral indices fixed. After an initial fitting, weak point sources {TS < 50) 
are removed. New point sources are identified by subtracting the best fit model from the 
counts map, then smoothing this gamma-ray residual map and identifying regions exceeding 
3 standard deviations above the average residual. 


Two additional components common to any model include the isotropic and Galactic 
inverse Gompton emission. The isotropic emission, originating from extragalactic diffuse 
gamma-ray emission and misclassihcation of cosmic rays in the LAT, is modeled by the 



Finally, we model diffuse emission arising from interstellar gas and dust. We use six 
different models for this analysis to determine the signihcance of the gamma-ray emission 
from each gas component. The baseline model, against which we compare all others, contains 
all the neutral gas elements: HI, GO, and dark gas {CODG). To check the signihcance 
of gamma-ray emission from the molecular cloud, we compare the baseline model to one 
containing only HI {HI). Two additional models are used to check the signihcance of gamma- 
ray emission from GO and dark gas individually: one model with HI and dark gas {DG) and 
another model with HI and GO (GO). To check whether the gamma-ray emission comes 
from an extended source or a point source, the GO and dark gas templates in the model 
are replaced with a point source located at the peak of the GO emission {PS). Finally, we 
test for any background point sources by adding one at the peak GO emission in addition to 
HI, GO, and dark gas {GODGPS). In the cases where the ROI has confused 2FGL sources 
(see table ED, we add the specihc 2FGL source into the model instead of adding a separate 
source. 


For neutral atomic gas, we separate the HI data from the LAB survey flKalberla et ah 


2 OO 5 I ) into two templates. In the hrst, we integrate over the velocity axis between ±20 kms 


-1 


to account for nearby, ambient gas associated with the cloud. The second template uses the 
rest of the HI data to account for background HI gas. We calculate the HI column density, 
A^(HI), using a spin temperature Ts = 125 K, consistent with previous gamma-ray studies 


^Sourse code can be found at https://sourceforge.net/projects/galprop 
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f Abdo et al. 2009 . 2010: Ackermann et al. 2012b ^0) 


JV(H I )(„. Ts) = -CTs In (l - . 


( 2 ) 


where C = 1.83 x 10^® cm“^ K“^, and T^g is the background temperature at 1.4 GHz. The 
A^(H 1) map is in a local coordinate system centered on the clouds and is linearly interpolated 
to 0?1 per pixel from the all sky map. 

For the H 2 template, we use IFco maps taken from Planck with a 2a cut. We use the 
“type 2” map, which combines data from different frequency channels to improve overall 
signal-to-noise by removing contaminating signals from CMB, dust, and free-free emission, 
though contaminati on remains. The contribution is compensat ed for by dividing 


kFco by a factor of 1.16 nPlanck Collaboration I l2014bl: IChen et al. 1120151) . 


In order to make a dark gas template, we remove a linear combination of iV(Hl) and 
Wc. Q from the visual extin c tion (A y) adapted from the map of color excess, E{B — V), from 
the \Planck Collaboration fl2014al) dust model, assuming Av=RyE{B — V), where we use 
the standard Ry = 3.1 here flCardelli et al. Ill989l: IWhittet et al. II 2 OOII) . This residual map. 


hereafter referred to as Ay res, is estimated as: 


Ay,res = ^ [n{R I ) -f /3WcO ) • 


( 3 ) 


The p arameters a and 13 are varied until the best fi t combination of A^(HI) and IFco is 


found fjCrenier et al. II2005I: lAckermann et al. Il2012al) . Both the Wqo and E{B — V) maps 


are use the same coordinates as the iV(Hl) map and are interpolated from the HEALPix 
maps provided by Planck. The four gas templates used for MBM 12 are shown in Figure HI 
These templates are all made 3° larger on each side than the region selected to account for 
possible gamma-ray photons leaking into the selected region due to the large point-spread 
function. 


Finally, the Fermi bubbles flSu et al. II 2 OIOI) represent an additional background source 
of gamma-ray emission for one cloud in this study, MB M 36. We add a constant component 
with a power law spectrum across the MBM 36 ROI fjSu et al. 1120lOl) to account for this 


emission. 


Every model component adds linearly to the total gamma-ray intensity flAckermann et ah 


2012al): 
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b) = Ej PS, + Ciso/iso + Ell i){e, b) 

+qco^co {(-lb) + gAv,res ^v,res (,^,b) + cic/ic(^5 b), (4) 

where Jic is the inverse Compton contribution to the gamma-ray intensity, Jiso is the isotropic 
contribution, and PS is the contribution due to point sources. The £t parameters are Qso and 
cic, which are normalization coefficients, and gHi,<?cO) and 5'Avres which are the emissivities 
of the respective gas templates in units of photon flux per H-atom cm“^, per K kms“^, and 
per magnitude of residual extinction, respectively. 

To perform the htting, we used the Python likelihood analysis tools provided by the 
FSSC. Using gtlike, we hnd an approximate solution with the ‘DRMNFB’ optimizer, and 
then rehne the solution by hxing point sources with TS < 100, removing any source with 
TS < 50, and rehtting with the ‘NewMinuit’ optimizer. 


4. Results 

We hnd diffuse gamma-ray emission from molecular gas to high signihcance in all 9 
regions studied. The control region has no diffuse gamma-ray emission from CO or dark 
gas. All clouds show extended gamma-ray emission, and all four 2FGL associations listed in 
Table [2] are identified as diffuse gamma-ray emission from the clouds themselves. Only MBM 
12 shows evidence of a background AGN. MBM 12, given its provisional classihcation as a 
blazar in the 2FGL and potential association with a radio source, goes through additional 
analyses discussed in a later subsection. 

Table [3] lists the TS values for each model compared to the baseline model CODG. Each 
model was £t separately to the data. The signihcance of gamma-ray emission from both GO 
and dark gas is given as TS'hs by comparing model CODG to HI. The signihcance of gamma- 
ray emission from the dark gas template is given as TS-dq by comparing model CODG to 
CO. The signihcance of gamma-rays from GO-emitting gas is given as TSqo by comparing 
model CODG to model DC. Determining whether any gamma-ray emission comes from an 
extended source is given as TSex by comparing model CODG to model PS. Finally, we test 
for contributions due to a background point source, presumably an AGN, given as TSagn, by 
comparing model CODG with model CODGPS. A US' > 20 indicates a significant diherence 
between the two models tested. 

In the control region, T^Ha = 1, indicating that we do not detect gamma-rays from 
gas traced by GO or dark gas. For the ROIs around the clouds MBM 04 and DIR 071-43 


there is no evidence for significant gamma-ray emission from gas traced by CO, as shown 
by a low TSqo- These ROIs show evidence for gamma-ray emission dominated by dark 
gas. Conversely, the ROIs for MBM 55 and MBM 32 show signihcant gamma-ray emission 
from gas traced by CO but not by dark gas; the TSqo is high while the T^dg is low. 
The remainings ROIs show evidence for emission from both gas traced by CO and dark 
gas to different degrees. For example, the ROI around MBM 12 show significantly more 
gamma-ray emission from gas traced by CO as opposed to dark gas, while the ROI around 
MBM 36 shows more gamma-ray emission from dark gas. Since this sample contains clouds 
dominated by gamma-ray emission from CO and some dominated by emission from dark 
gas, the translucent clouds may represent the transition to fully molecular clouds. 

For every ROI, TRex is large enough to conclude that the gamma-ray emission is not 
coming from a single point source. A TSagn > 20 indicates the presence of a background 
AGN. Only MBM 12, discussed in the following section, shows significant evidence for a 
background source. Thus, out of the four point sources given in Table [21 we confirm the 
existence of one background gamma-ray point source and attribute the other three to diffuse 
gamma-ray emission from molecular clouds. 


We also compare the baseline model to one utilizing the Fermi standard G alactic diffuse 
model for Pass Tjj. However, only one cloud in our survey had GO data from iPame et ai 


0200 ll) included in the diffuse model (MBM 12). Therefore our models, which include GO, 
should more accurately reflect the gas origin of the gamma-ray emission. Further, as will be 
discussed in Section 14.21 the dust alone results in a poor £t to the data. 


Over the entire energy range, when our gas templates all assume a single power law 
spectrum, the standard diffuse model has a higher likelihood than any of the models we 
fit {TS > 400). The higher likelihood of the Fermi Galactic diffuse model is a result of 
it being fit in many small energy bins, rather than assuming a single power law over the 
entire energy range, as we do. A single power law does not capture the curvature in the 
spectrum. Therefore, we tested a log-parabola and a broken power law spectrum for the HI 
template and find the fits strongly favored the broken power law spectrum with aTS of at 
least 200 over the single power law. The Fermi Galactic diffuse model still £t better with 
Si TS ^ 100. However, for MBM 12, we have enough data to fit several energy bins, more 
closely replicating the Fermi Galactic diffuse model. In every energy bin, our model fits the 
data as well or better than the standard Galactic diffuse model. Thus, we validate our model 
based on the Planck GO, LAB HI, and Planck dust data. 


For the broken power law spectrum, the average energy break occurs at 1.19±0.10 GeV 


^gaL2yearp7v6_source 





9 


and the power law index changes from —1.80 ± 0.16 to —2.82 ± 0.06. The break at around 1 
GeV is expected from cosmic ray proton interactions with the cloud and the power law index 
above the break energy agrees with predict ions for garnma-ra y emission from a cloud located 
far from any cosmic ray acceleration site flAharonian 120041) . Changing the spectral shape 
from a single power law to broken power law results in a 30% increase in the HI emissivity 
but does not affect the other gas emissivities. 


We leave the CO and spectra as single power laws for simplicity because HI 

is the dominant source of gamma-rays across the ROI. Changing the spectrum of the HI 
CO, and Hy^res templates to a broken power law resulted in no improvement (TS < 10 for 
six additional free parameters) over the model where only the HI had a broken power law 
spectrum. 


4.1. MBM 12 


MBM 12 is the best studied high latitude, translucent molecular cloud fe.g.. 

Pound et ah 

1990; 

Ingalls et ah 

1994; 

Moriartv-Scieven et ah 1997: 

Timmermann et ah 

19981. It has 


the highest peak ITco of the sample and was tentatively identihed as a blazar with a “con¬ 
fused” designation in the Fermi 2-year catalog (2FGL J0257.9-|-2025c). It is also coincident 
with a radio source from the 3'’'^ MIT-Green Bank radio survey (MG3 J025805-I-2029). MBM 
12 , therefore, merits further analysis. 

Figure 0] shows the distribution of the gas in the region around MBM 12. Figure [5] 
shows the total number of gamma-rays detected in the region around MBM 12 over the time 
period analyzed. A 2FGL source, as well as the coincident radio source, he near the centroid 
of the molecular cloud. 

To establish the molecular cloud origin of the gamma-rays over an AGN source, we 
perform both a spatial and a variability test. The results of the spatial tests are given as 
TSag'n and TSe^- MBM 12, after removing 2FGL J0257.9-|-2025c, has TS'agn = 96 which 
verihes the presence of the 2FGL point source. All other T S values are calculated with this 
point source included in the model. A TS'ex = 308 implies that the gamma-ray emission is 
extended. The results from these two spatial tests argue in favor of the translucent cloud 
origin of the gamma-ray emission. 

For the variability test, we compare the gamma-ray flux in 15-week time bins to the flux 
over the entire 58 month period of the survey. As MBM 12 lies far from possible cosmic ray 
acceleration sites, the cosmic ray source is dominated by the mean steady-state diffusion of 
cosmic rays through the Galaxy, and should be constant. As a result, signihcant variability 
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in the gamma-ray flux from MBM 12 would indi cate a background AGN, which are the 
most comm on gamma-ray sources at high latitude flNolau et ah II2012I) and typically exhibit 


variability flAckermaun et ah Il2011bl) . 


The lightcurve of MBM 12 is given in Figure [6] and is the sum of the gamma-ray 
fluxes from HI, CO, and dark gas. By eye, they all appear to have constant emission over 
the observed period. Time bins 4 and 5 as well as 10 and 11 were combined to minimize 
convergence errors. The dashed line corresponds to the integrated flux found over the entire 
time period, and the shaded region is the associated uncertainty. With 13 time bins, we And 
Xi 3 = 1.1 when comparing the flux to a constant, implying that the gamma-ray flux from 
MBM 12 is consistent with zero variability. 


The tests show that we both detect diffuse gamma-ray emission from MBM 12 and 
confirm the existence of 2FGL J0257.9-|-2025c, associated with MGS J025805-I-2029. Because 
we can recover point sources that have been intentionally removed from the model, we gain 
confidence that we can detect other new sources not included in the 2FGL and disentangle 
them from diffuse cloud emission. 


As the brightest source in the survey, we also use MBM 12 to test the validity of our 
chosen energy range. Data analyzed in this region between 10 to 100 GeV shows no significant 
emission from molecular gas, with a T = 4. Additionally, the HI emissivity between 250 
MeV and 100 GeV differs by only 5% from the HI emissivity between 250 MeV and 10 GeV. 
This difference is well within the systematic uncertainty. Extending the analysis to 100 GeV, 
therefore does not significantly change the emissivity. 


4.2. Uncertainties 


We are ultimately trying to probe for small effects, such as possible gradients in the 
cosmic ray flux in the Solar neighborhood, so it is important to identify and characterize 
as many sources of uncertainty as possible. Systematic errors introduced by the LAT in¬ 
strument are estimated to be around 10% (INolan et ah Il2012l) . They arise primarily from 
uncertainties in the instrument response function, the energy determination of the photons, 
and the effective area of the LAT. The remaining uncertainties are derived from the results 
of the likelihood analysis. 


Non-local HI, GO, and dark gas each have two free parameters associated with them 
and local HI has four free parameters. In addition, the inverse Gompton and isotropic 
emission templates include a normalization term, bringing the total number of template free 
parameters to 12. Additional, compounded, free parameters are used to create the Ayres 
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maps and are discussed below. Point sources add additional free parameters and models 
with a lot of free parameters tend to converge poorly or force a parameter to one of its 
limits, which skews the error calculation. Therefore, our analysis procedure removes weak 
point sources. We test models with different initial values and find that when the models 
converge, the resulting gamma-ray fluxes and photon indices typically vary by less than 10%. 
We take the combined uncertainty due to the LAT systematics and the likelihood variations 
to be ~15%. 


The uncertainties in the emissivities depend on the detailed model inputs. The LAB 
survey measured the radiatio n roughly —400 kms~ ^ to -1-400 kms“^ around the 21 cm line 


with a sensitivity of 0.09 K flKalberla et ah 11200511 . The high velocity gas is all very far 


away while the low velocity gas is much closer, so we perform a velocity cut of ±20 kms“^on 
the data to separate nearby gas from far. Most of the HI emission in our clouds lies inside 
of this cutoff. The velocity cutoff is a fairly small source of error. Generally, the higher 
velocity gas contributes at most 10% more column density. High velocity gas, given its much 
farther distance, is expected to contribute little to the total gamma-ray flux. When the 
higher velocity HI is significantly detected, the effect on the HI flux of the cloud is < 10%. 

In addition, the HI template suffers from uncertainty due to the assumed spin temper- 
ature. Spin tempera tures likely change across an ROI and even across an individual cloud 


flFukui et ah Il2014al) . We tested four spin temperatures: Ts = 80 K, 125 K, 400 K, and 


Ts —)■ cxo. The gamma-ray flux from HI changed less than 7% while the emissivity decreased 
by 15% with increasing Ts over the entire range of Ts values. 

Surrounding the molecular cloud may be a shell of optically thick HI with a very low 


spin temperature, that a lso contributes to the dark gas phenomenon fiFukui et ah 12015 


Stanimirovic et ah 1120141) . Any unaccounted for HI due to a lower spin temperature than 
the 125 K used should be captured in the Ay,res rnap. This should not affect the HI emissivity 
or gamma-ray flux by more than 5%, as the HI in the cloud accounts for less than 5% of 
the total HI emission across the ROI. 


The main uncertainties from the CO template comes fro m the 2cr cutoff used to remov e 
noise, and the ^^CO contamination of the Planck IFco data ^Planck Collaboration 1 12014bl) . 
We tested a Icr cutoff as well and saw emissivities and fluxes systematically lower by around 
5%. This is expected since the CO covers a larger area in the lex template versus the 2 ct 
template. Yet the distribution of CO does not change much; most of the difference between 
the Icr and 2cr templates is noise and distributed roughly uniformly around the map. Nearly 
the same flux is being emitted from a larger amount of CO, therefore the emissivity decreases 
accordingly. 
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The Avrpfi map is a li n ear com bination of the LAB HI map, the Planck CO map, and 


the \Planck Collaboration fl2014aj) color excess map, where the nncertainties for the color 
excess are generally less than 8 %, with an average across each ROI of less than 4%. The 
nncertainty in the color excess dominates the Ay,res errors over those of the HI and CO 
maps. The fluxes and emissivities for CO and Ay^res are insensitive to the change in spin 
temperature. Even between the two extremes, T 5 = 80 K and T 5 —)■ 00, the Ay^res template 
changes by less than 2 % on average, which is much smaller than the uncertainties of the 
color excess map. 

Additional uncertainties in the Ay^res rnap arise in Ry = Av/E{B — V). While overall 
variations in Ry will not affect the HI emissivity, CO and dar k gas emissivities ma y be 


affected if Ry varies across the ROI. Ry = 3.1 in diffuse regions fICardelli et ah 


increases to Ry > 4 in some molecular clouds fjVrba et al 


1993 


Kandori et al. 


■9891) but 

2 nnl . A 


constant Ry may underpredict the extinction values in molecular clouds by as much as 30%. 
The uncertaintiy in the emissivity is not clear due to the relationship between Ay^res and 
IRco- 

To study this effect, other dust templates may be more appro prate. For example, a 


proper ly scaled T 353 map is recommended for regions of higher Ay UPlanck Collaboration 


2014al) instead of the E(B — Rlmap. This test was performed in the Chamaeleon region in 
Fermi /Planck Collaboration fj2014d ). Changing the dust map necessarily affects the Ay^res 
map, and therefore also the htted and subsequently calculated quantities. We test both dust 
maps in MBM 12. Comparing both models yields a.TS = 2. Therefore neither E{B — V) 
nor r 353 are preferred over the other. The Ay^res emissivity changed by up to 30% while 
the HI and CO emissivities changed by less than 5%. We quantify the effect of the dust 
template more thoroughly in the full survey and defer that discussion to future work. We 
argue that with only 9 clouds in our current sample, a 30% change in dark gas emissivity 
does not signihcantly affect the results discussed in Section El 

Finally, we also tested a model using the dust map as the sole gas tracer. In prin¬ 
ciple, the extinction should trace all gas species and thus should provide a model with 
signihcantly fewer free parameters. The dust model, however, does not reproduce the 
gamma-ray data. The TS between model CODG and the model with dust alone is 175; 
the dust model hts signihcantly worse than our baseline model, and exhibits large, struc- 
tured residuals in regions with CO ari d dark gas. This may be due to varying dust properties 
Fermi / Planck Collaboration Il2014d) . a lack of sensitivity in dense regions, a limited range 
of applicablility, or a nonlinear response. This validates our combination of gas templates. 


In sum, there a is 15% systematic uncertainty in the gamma-ray hux from both the 
LAT instrument and the likelihood analysis along with an additional uncertainties in the 
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emissivities due to the choice of model of , Qco ± 5%. For this study, we adopt 

gAv.res ± 8%. 


Discussion 


From these observations, we may derive the gamma-ray emissivity, the gamma-ray pho¬ 
ton flux per H-atom, for the brightest clouds within 270 pc of the Sun between 250 MeV 
and 10 GeV. The £t treats the energy range as a single bin. Table H] lists the emissivities 
for HI, CO, and the dark gas, while Table E] lists the remaining model parameters exclud¬ 
ing point source normalizations. The HI emissivities are all generally consistent with each 
other, and close to the average emissivity of gHi(250 MeV — 10 GeV) = (8.1 ± 1.4)xl0“^'^ 
photons cm“^ s“^ sr“^ H-atom“^. This result is consistent with the emiss ivities, all between 

n the Cygnus region, (8.76 ±0.33) x 10“^^ (j Acker maun et ah 


2012a 

), and that found in 

Ackermann et al. 

f2012c) 


1-27 


()2012d) . The Chamaeleon region has be en analyzed twice and th e HI 


(10.8 ± 0.4) X 10-2R 

Fermi /Planck Collaboration 1 

2014 

c) and the HI emissivity around R 

Coronae Australis was found to be (10.2 ±0.4) x 10 

-27 

Ackermann et ah 

2012c 

). All values 


are less than 2 a from our results. 

The emissivities for the molecular gas tracers, the CO and dark gas, vary more sig- 
nihcantly. The CO component has an average emissivity of geo (250 MeV - 10 GeV) = 
(1.6 ± 0.6)xl0“® photons cm“^ s“^ sr“^ (K kms“^)“^ and the dark gas component has 
an average emissivity of gAvres(250 MeV - 10 GeV) = (1.5 ± 0.7) x 10“® photons cm“^ 
s“^ sr“^ mag“^, taking into account statistical errors. The CO emissivity was found to be 
(3.01 ±0.16) X 10“® in the Cygnus region while the clouds Chamaeleon, R Coronae Australis, 
and Cepheus/Polaris flare have (1.04±0.08) x 10“®, (1 .9±0.2) x 10~^, and (1.23± 0.05) x 10“® 


photons cm ^ s ^ sr ^ (K kms ^) respectively llAckermann et ah l2012crl. The dark 


gas emissivity for all four regions was (2.75 ± 0.26) x 10 ^ flAckermann et ah Il2012al) and 
(1.36 ± 0.04) X 10-^ (2.3 ± 0.2) x 10-^ and (1.38 ± 0.08) x lO"® photons cm'^ 


sr 


mag ^ (jAckermann et ah ll2012cD . respectively. Our results are entirely consistent with the 


previously found values. 
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5.1. Cosmic Rays in the Solar Neighborhood 

A change in the gamma-ray emissivity indicates a change in the incident cosmic ray flnx. 
On the right of Fignre[9l we plot the HI emissivity and see no variation. MBM 12 has a high 
valne, thongh systematic nncertainties place it within the 2 a range of th e average HI emis- 
sivity. A constant HI emissivity is in agreement with the conclusions o f I Acker man n et ah 


fl2ninli 


(]2012cll and the predictions of GALPROP models fIStrong et ah 1120041) . lAbdo et ah 
looked for gradients in gni and compared it to predictions from GALPROP. They measure 
a 10% decrease in gni from the Gould Belt to the local arm, a distance of about 1 kpc. The 
maximum Galactocentric distance between any of the clouds in this survey covers about 0.25 
kpc, so we might expect to see a 2.5% change in gni across our sample, which is well within 
our quoted uncertainties. 

Additionally, we report the gamma-ray photon index from the GO to isolate the molec¬ 
ular cloud. The indices from GO and dark gas components are identical within the uncer¬ 
tainties. This index is similar to the average of the two power law indices in the broken 
power law spectrum of HI. The indices of the clouds are listed in table H] and are taken 
from the GO fit. Where GO is not detected, the index is taken from the dark gas fit. The 
left side of Figure [9] shows the power law index for every cloud. They all he close to the 
best fit value —2.25 ± 0.10 with no evidence of a variation. The lack of detected variations 
implies a uniform cosmic ray flux incident on each cloud verifying that there are no cosmic 
ray sources near the clouds studied. The farthest a cloud can be from a supernova remnant 
and still recei ve a cosmic ray excess is 100-200 pc for a 10^ year old supernova remnant 


flGabici II 2 OIII) . 


5.2. X-factors 


We trace molecular gas with two components, Wco and Ayres, and so we require two 
conversion factors to estimate the column de nsity of H^ . Ayres is not entirely molecular gas, 


however; sorne fra ction is atomic hydrogen flFukui et ah Il2014al: IStanimirovic et ah II2014J: 


Fukui et al. II2015I) . The molecular gas column density can therefore be written as: 


7V(H2) = ^CO + /Xi^Ayres , (5) 

where / is the molecular fraction of the dark gas. This expression is proportional to the 
traditional X-factor, Xco, where GO is assumed to trace all the H 2 . 


Analysis of gamma-rays alone does not determine what fraction of dark gas is H 2 , but 
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we can put upper and lower limits on Xco: 


X, 


CO 


= 


I f^V,Tes \ 

V w^co ; ■ 


( 6 ) 


Assuming a constant cosmic ray flux and that cosmic rays penetrate the entire cloud, 
every proton should be subject to the same num ber of cosmic ray interactions. This as sump- 
tion was verifled in the Chamaeleon clouds in \ Fermi / Planck Collaboration fl2014c|) . and 
implies gH 2 = 2gHi- The gamma-ray emission is proportional to the number of molecules, 
but we use Wco in our analysis: gH 2 X(H 2 ) = qco^co ■ This leads to the relationship 
Xqq = gco/2gHi- A similar argument leads to the relationship X^^ = (?Avres/([l + /I'lni), 
where the factor [l + f] arises because a fraction of the gas represented by molecular hydrogen 
which has two protons. 


From the emissivities of MBM 12 in different energy bins given in Table [6l we plot qco 
versus gni in each energy bin in Figure [71 We And a linear relationship between the two 
emissivities and plot the best At line. Similarly, we plot qAvres versus gni in Figure [HI Again, 
we see a linear relationship, which gives us confidence that the cosmic ray flux at MBM 12 
is constant and penetrates through the entire cloud. For MBM 12, the slope of the best fit 
line in Figure [Tlis X^q = (4.8 ± 1.2)xl0^® cm“^ (K kms“^)“^, which agrees with the value 
in Table [7[ obtained as a ratio of the emissivities from Table [H The slope of the best fit line 
in Figure [H] is Xj^^ = (22.4±3.5)x 10^° cm“^ mag“^, roughly 2.5 ct from the value in Tabled 

Most clouds are not detected significantly enough to make a spectrum, so we cannot 
verify the linear relationship in each case. MBM 12 has the highest Ay of any of the clouds 
in this study, so we conclude cosmic rays penetrate through every cloud in this study. That 
our clouds are far from cosmic ray acceleration sources and the HI emissivity is constant 
across them gives us some confidence that the cosmic ray flux is constant across an entire 
molecular cloud. We therefore assume X^q = gco/2(3'Hi and = 5'Avres/([l + /I'Ihi) for 
every cloud in this study. 


The average value for X^q among the clouds is (1.1 ± 0.4)x 10^° cm ^ (K kms ^) ^ 


This is consistent with previous gamma-ray studies of nearby molecular clouds ([Ackerman n et ah 


2012cl) . where Xco ^ 1x10^° cm“^ (K kms“^)“^. This value is lower than that found 
for high latitude clouds of 1.67x10^° cm“^ (K kms“^)“^ ( Paradis et ah 2012 1. Our X'qq 
for Cham-East II is consistent with a recent analysis of the Chamaeleon cloud complex 
\FermilPlanck Collaboration [[2014c[) which finds Xco ~ 7x10^® K kms“^. However, X'^^ 
is always higher than X^q as seen in Table [TJ with an average of (19.0 ± 8.8) x 10^° . As long 
as any fraction of Ay,res represents H 2 , the combination of the two will increase Xco- 
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Assuming / = 1, we report the average Xco of each cloud in Table [71 The average among 


-2 


(K kms ^ with large c 


oud-to-cloud varia- 


all the clouds is Xco = (1.6 ± 0.5) x 10^° cm 
tions. This result agrees w ith previous studies of Xqq at high latitudes fjMagnani fc Onello 


19951: [Paradis et ah II2012I) . suggesting the dark gas in these previous studies i s primarily 


mole cular. Our average Xqo is higher than the average value found in Perseus flLee et ah 


2014[ ). which may suggest that the molecular fraction of the dark gas in Perseus is signif¬ 
icantly lower than 1 or cosmic rays do not penetrate deeply into the CO-bright regions of 
the cloud. Our Xco for every cloud is consistent with that found in th e high latitude cloud 
MBM 40 which found an average Xco = 1.3x10^° cm“^ (K kms“^)“^ f Gotten fc Magnani 

2018h . 


Figure [10] shows Xco as a function of Galactocentric distance. The solid line in the hgure 
shows one of the more extreme variations of Xco considered near the Solar neighborhood 
(llsrael [[2000h . As with the gamma-ray spectrum, we find no evidence for overall variation 
in Xco over this small extent of 270 pc. 

Figure [m shows Xco across the CO-bright part of MBM 12, assuming / = 1 uniformly 
across the entire cloud. This places an upper limit on the value of Xcq. 77co appears to 
increase towards the edge of the CO-emitting region where the total extinction drops by as 
much as an order of magnitude. The middle of the cloud shows a low Xco, where the total 
extinction increases due to increased gas density. A higher density supports the transition to 
fully molecular gas. This transition reduces the fraction of dark gas in the cloud. Thus, the 
departure of Xco from X^q in the CO-bright part of the cloud may reflect the transition of 
the cloud from atomic to molecular gas or it may reflect a level of small-scale dumpiness in 
the cloud, supporting some interior dissociation along the line of sight. If / decreases to zero 
toward the edge of the CO-bright part of the cloud, we can potentially r ecover a constant 
Xco across the cloud, though not expected to occur fjWolfire et ah [[2010[) . 


As can be seen by visual inspection of the CO and dark gas maps i n Figure [H dark gas 


extends beyond the CO boundaries as expected from numerical studies flWolhre et ah [[2010 


Velusamv et ah [[2010[) . In these regions, we can only place limits on Xcq. Therefore, we do 
not attempt to estimate Xco beyond the CO-bright boundary of the molecular cloud. 

Finally, as an initial evaluation of the dust model, we replace the Planck E{B—V) map in 
MBM 12 wi th the Planck dust opacity map, scaled by the given ratio: E{B — V)/t 353 = 
1.49 X 10^ f Planck Collaboration [2014al l. In MBM 12, changing the dust model greatly 
affects X^^. However Xco changes by ±0.5x10^° cm“^ (K kms“^)“^, which is smaller than 
the variations of Xco. is thus moderately insensitive to large variations in X)^^ due to 
suppression by the ratio Ay,res /hFco , which averages between 0.01 - 0.10 mag (K kms“^)“^. 
It is worth noting that, while the distribution of Xco seen in Figure [TTl depends on the dust 
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tracer used, the magnitude of the variation of Xco is larger than any uncertainty of Xco- 
We will quantify the uncertainty due to changing the dust tracer in future work. 


6. Conclusions 


We study the gamma-ray emission from nine high latitude, translucent molecular clouds 
and hnd the gamma-ray spectrum and emissivity from the gas in the clouds. All nine 
Planck-selected CO clouds were signihcantly detected, showing extended emission associated 
with molecular gas. We estimate the systematic uncertainties associated with modeling the 
gamma-ray emission with HI, CO, and dark gas templates. The HI emissivity does not vary 
across the regions in the sample. For some clouds, the gamma-ray emission is dominated 
by the CO-emitting gas, while for other clouds the CO-dark gas dominates the gamma-ray 
emission. 


Xqq ranges from 0.2x10^° cm ^ (K kms ^ to 1.6x10^° cm ^ (K kms ^ with 


21 2 1 order 


large uncertainties, and spans a similarly large range around 10^^ cm“^ mag 
to compare these X-factors to the traditional conversion factor between Hco and H 2 Xqq 
and are added together assuming the molecular fraction of the dark gas is 1. The results 
are consistent with previous studies and suggest no change in cosmic ray flux across ~ 300 
pc region around the Solar Syst em. The combinat i on of and X^^ may also explain 
the low values of Xno found in lAckermann et ah (j2012d ) compared with estimates using 
other methods flBolatto et ah 1120131) . However, the choice of dark gas tracer should be made 
carefully. 
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£= 180 ° ^— X (pc) — >e=o 


Fig. 2.— Same as Figured! except the vertical direction indicates distance from the galactic 
plane and the Galactic Center lies off to the right. 
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£=270° ^— X (pc) —>£=90° 


Fig. 3.— Same as Figured! except the vertical direction indicates distance from the galactic 
plane. 
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Fig. 4.— Gas templates for MBM 12. top-left. Local HI column density in 10^^ cm“^, top- 
right. non-local HI column density in 10^° cm“^, bottom-left. IFco in K km s“^, bottom-right. 
the dark gas template, Hyres in magnitudes. 
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Fig. 5.— MBM 12 gamma-ray counts map, smoothed with a a = 0.5° Gaussian, with CO 
contours at Wco = 5,10,15,20, and 25 K kms“^. The cyan circles mark the positions of 
signihcant 2FGL point sources, and the black x marks the position of a radio source behind 
MBM 12. The counts map covers a smaller area than the gas templates in Figure 01 
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Week Number (MET) 


Fig. 6.— Lightcurve of the total gas flux from the MBM 12 ROI in mission elapsed time 
(MET). The dashed line is the flux determined by htting the ROI over the entire time range 
and the gray shaded region is the statistical uncertainty on this flux. 






Fig. 8.— Gamma-ray 
linear relationship. 
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H-atom“^) 


for dark gas versus emissivity for HI in MBM 12 shows a 
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Galactocentric Distance (kpc) 


Fig. 9.— Power law index for gamma-ray emission from CO (left) and HI emissivity (right) 
as a functi on of Galactocentric distance, with the distance to the Galactic Center Ro = 
8.33 ±0.35 flGillessen et al 1120091 ). marked with a vertical dotted line. The horizontal dashed 
lines are the best £t constant power law index and HI emissivity, respectively, and their 
uncertainty is the shaded region. 
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Fig. 10.— Xco as a function of the distance from the Galactic Center, as in Figure [9l 
The solid lin e represents a variati on in Xco due to metallicity from llsrael (120001) . rescaled 


according to IPineda et ah (120131) 
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Fig. 11.— A map of Xqq from Equation [6] assuming / = 1, across the CO-emitting portion 
of MBM 12 in units of 10^° cm“^ (K kms“^)“^. 
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Table 1. Cloud Properties 


Name 

1 

(deg) 

b 

(deg) 

Peak lEco 
(K kms“^) 

Distance 

(pc) 

Area 

(deg2) 

MBM 36’" 

3.9 

35.7 

27.85 ± 0.74 

105 ±7 

4.041';^® 

DIR 071-43 

71.6 

-42.8 

3.69 ±0.70 

170 ± 20'" 

1 Qc^+i.e 
-‘-•'^""-0.85 

MBM 55"= 

88.5 

-41.2 

15.75 ±0.70 

2061^ 

21.23li®7f 

MBM 02 

108.3 

-51.9 

7.45 ± 0.73 

20Qt]^ 

1 '77+0.78 
-*-• ' ' -1.64 

MBM 04^ 

132.9 

-45.6 

15.31 ±0.78 

2Q9tll 

1 1C+9.6 

-‘-•■‘-'9-0.39 

MBM 32*^ 

147.6 

40.7 

8.06 ±0.62 

269^« 

2.36lli 

MBM 12f 

159.1 

-34.3 

29.78 ± 0.73 

234i“ 

io.24i«:?j 

MBM 20® 

210.9 

-36.5 

12.39 ±0.72 

1241“ 

2.48lJ:«t 

Cham-East iC 

312.9 

-28.6 

20.25 ±0.53 

150 ± 20'" 

12.45+°^ 

Control 

250.0 

30.0 

< 2.46 ±0.7 

— 

< 0.06 


Note. — Regions around the Wco peak included in the gamma-ray 
analysis includ e addi tional molecular clouds. Distances derived from 


Schlaflv et al. (1201411 unless otherwise noted. 


'^ROI also includes: MBM 33 through MBM 35, MBM 37, and MBM 
38 


b 

c 


Distances derived from iLallement et al. (1201411 
ROI also includes: MBM 53 and MBM 54 


‘^ROI also includes: MBM 03, DIR 121-45 

®ROI also includes: MBM 27 through MBM 31, HSVMT 24, 
HSVMT 27, and HSVMT 28 

^ROI also includes: MBM 07 through MBM 09, and MBM 11 
®ROI also includes: DIR 203-32 


^ROI also includes: Cham-East I 
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Table 2. Coincident Point Sonrces 


Cloud Name 

Fermi Oatalog 
Name 

Power Law 

Index 

Variability 

Index®^ 

Olassihcation 

Flags 

(AT) 

MBM 12 

2FOL J0257.9+2025C 

-2.19 

25.55 

blazar’=’ 

6 

MBM 36 

2FOL J1553.5-0324 

-2.22 

15.85 

— 

— 

0313.1-28.6'= 

2FOL J1942.7-8049C 

-2.40 

32.5 

— 

2,6 

0315.1-29.0'= 

2FOL J1925.7-7836C 

-2.31 

14.1 

— 

1,6 


Note. — The molecnlar clonds coincident with point sonrces in the 2FGL. Properties 
are from the 2FGL. A bla nk means there is n o classihcation or flag dnring the htting 
process. Relevant flags are flNolan et ah 1120121) : 

= 1: Source not detected signihcantly when the diffuse model was changed. 

N = 2: Source location changed beyond its 95% error ellipse when the diffuse model was 
changed. 

N = 6: On top of an interstellar gas clump in the model of diffuse emission. 

^'A source with a variability index above 41.6 is considered variable at the 99% conh- 
dence level ( Nolan et ah 2012 1. 


Nolan et ah (1201211 designation: bzb, tentative blazar classihcation. Coincident with 


the radio source MG3 J025805+2029 


Part of Cham-East II 















Table 3. 


Model Test Statistic Values 


ROI name 

TSn, 

TSdg 

TSco 

TRex 

TS'agn 

MBM 36 

595 

140 

92 

314 

9 

DIR 071-43 

76 

47 

2 

76 

7 

MBM 55 

341 

0 

273 

318 

-2 

MBM 02 

62 

10 

17 

41 

1 

MBM 04 

80 

47 

18 

80 

18 

MBM 32 

166 

9 

76 

122 

6 

MBM 12 

691 

37 

346 

308 

96 

MBM 20 

457 

43 

63 

124 

0 

Cham-East II 

832 

93 

453 

624 

16 

ControR 

1 

0 

-1 

- 

- 


Note. — The significance of various diffuse emission 
components, rounded down to the nearest integer. 

^In the control region, the models with point sources 
were not fit because there was no significant CO emission 
and little dark gas in the ROI. 
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Table 4. Cloud Parameters 


ROI Name qu i qco QAvrea Photon Index 


MBM 36 
DIR 071-43 
MBM 55 
MBM 02 
MBM 04 
MBM 32 
MBM 12 
MBM 20 
Cham-East II 
Control 


9.1±0.5l?:g 

5.9±0.3lJ;;^ 

8.9±0.9l2;0 

6.9±0.6lJ;® 
7.9±0.5lJJ 
7.1±0Atli 

12.4 ±0.2+1® 
7.7±0.4lJ;^ 
6.9±0.3li:^ 
9.0±0.5l^;^ 


1.2 ±0.1 ±0.2 

1.1±0.7±0.2 

1.4 ±0.1 ±0.2 

2.4 ±0.7 ±0.4 

2.1 ±0.2 ±0.3 

1.4 ±0.1 ±0.2 
2.8 ±0.4 ±0.4 
1.3 ±0.1 ±0.2 


3.2 

± 

0.3 

± 

0.5 

1.6 

± 

0.3 

± 

0.3 

0.3 

± 

0.2 

± 

0.1 

o 

bo 

± 

0.3 

± 

0.2 

1.9 

± 

0.4 

± 

0.3 

0.5 

± 

0.2 

± 

0.1 

1.3 

± 

0.2 

± 

0.2 

1.2 

± 

0.2 

± 

0.2 

4.3 

± 

0.5 

± 

0.4 


-2.13 ±0.10 
-2.17 ±0.44 
-2.23 ±0.06 
-2.40 ±0.19 
-2.34 ±0.IP 
-2.27 ±0.07 
-2.05 ±0.08 
-2.47 ±0.10 
-2.32 ±0.06 
-2.61 ±0.05'" 


Note. — Gamma-ray emissivity of gas templates for 250 MeV < E < 10 GeV 
with the associated statistical and systematic uncertainties. The photon index 
reported is from the CO template. 

Units: photons s“^ sr”'^ H-atom”'^), gco(10“® photons cm“^ s“^ sr”'^ (K 

kms“^)“^), gAvres(10~^ photons cm“^ s“'^ sr“^ mag“^). 

'"CO not detected signihcantly, index taken from dark gas. 

'"Neither CO nor dark gas detected signihcantly, index taken from HI above ~1 
GeV. 







Table 5. Other Fitted Parameters 


ROI Name 

gniJar 

CiC 

^iso 

MBM 36 

— 

1.2 ±0.7 

0.7±0.1 

DIR 071-43 

— 

2.6 ±0.3 

0.9 ±0.1 

MBM 55 

11.2 ± 1.8 

— 

1.3 ±0.1 

MBM 02 

— 

2.0 ±0.6 

1.1 ±0.1 

MBM 04 

— 

2.0 ±0.5 

1.0 ±0.1 

MBM 32 

12.4 ±1.2 

1.2 ±0.2 

0.8 ±0.1 

MBM 12 

— 

— 

0.7±0.1 

MBM 20 

— 

1.5 ±0.4 

1.0 ±0.1 

Cham-East II 

— 

— 

1.3 ±0.1 

Control 

13.6 ± 1.6 

— 

1.0 ±0.1 


Note. — Gamma-ray emissivities for the non¬ 
local HI and the normalization factors for the in¬ 
verse Compton and isotropic components for 250 
MeV < E < 10 GeV with the associated statistical 
uncertainties. Blanks indicate the component was 
not detected signihcantly and thus not included in 
the model. 

Units: gm f ar (10~^’^ photons s“^sr“^ H-atom“^) 







Table 6. MBM 12 Parameters 


Energies 

(MeV) 

gni 

Qco 

^-^V,res 

^iso 

250 - 400 

3.24 ± OAAtolt 

0.24 ±0.09 ±0.04 

1.46 ±0.26 ±0.25 

1.4 ±0.2 

400 - 630 

2.45±0.19^[;i^ 

0.30 ±0.05 ±0.05 

0.45 ±0.14 ±0.08 

1.3 ±0.2 

630 - 1000 

1.54±0.15l°i5 

0.21 ±0.03 ±0.03 

0.33 ±0.08 ±0.06 

1.3 ±0.3 

1000 - 1580 

0.88 ± 0.04+°;20 

0.08 ±0.02 ±0.01 

0.22 ±0.05 ±0.04 

1.0 ±0.1 

1580 - 2510 

0.52 ± 0.07+[j;J^ 

0.07 ±0.01 ±0.01 

0.09 ±0.03 ±0.02 

— 

2510 - 3980 

0.19±0.03l[;;[^^ 

0.04 ±0.01 ±0.006 

0.02 ±0.01 ±0.003 

1.4 ±0.4 

3980 - 10000 

0.08 ± 0.03l[;;[^? 

0.03 ±0.01 ±0.005 

0.03 ±0.01 ±0.005 

1.8 ±0.5 


Note. — Gamma-ray emissivities and the associated statistical and systematic nn- 
certainties for the gas templates. The isotropic component only inclndes statistical 
nncertainties. Non-local HI and inverse Compton are not detected signihcantly in any 
individual energy bin and are not included. 







Table 7. 


Calculated Cloud Properties 


ROI Name 

(cm“^ (K kms“^)“^) 

(cm“^ mag“^) 

Vco/lO^o 

(cm“^ (K kms“^)“^) 

MBM 36 

6.8 ±2.0 

17.6 ±4.4 

1.4 ±0.5 

DIR 071-43 

9.3 ±6.7 

13.6 ±4.0 

1.9 ±0.9 

MBM 55 

7.6 ±2.3 

1.7 ± 1.2 

0.8 ±0.2 

MBM 02 

17.5 ±7.0 

5.8 ±2.6 

2.3 ±0.8 

MBM 04"^ 

— 

11.8 ±3.8 

— 

MBM 32 

14.8 ±4.5 

3.5 ±1.6 

1.5 ±0.4 

MBM 12 

5.5 ±1.8 

5.4 ± 1.5 

0.7 ±0.2 

MBM 20 

18.0 ±5.6 

7.8 ±2.2 

2.5 ±0.8 

Cham-East II 

1.1 ±2.6 

35.3 ±9.1 

1.7 ±0.7 


Note. — Values for the conversion factors between Wqq and dark gas. The 
last column is the average value across the cloud of the combination of X'qq and 
X'j^^ as in Equation |5] assuming / = 1. Uncertainties include systematics from 
the emissivities. 


'^CO not detected signihcantly. 






